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The BEM with graded meshes for the electric held 
integral equation on polyhedral surfaces 

A. Bespalov* S. Nicaise^ 


Abstract 

We consider the variational formulation of the electric field integral equation on 
a Lipschitz polyhedral surface F. We study the Galerkin boundary element discreti¬ 
sations based on the lowest-order Raviart-Thomas surface elements on a sequence of 
anisotropic meshes algebraically graded towards the edges of F. We establish quasi- 
optimal convergence of Galerkin solutions under a mild restriction on the strength 
of grading. The key ingredient of our convergence analysis are new componentwise 
stability properties of the Raviart-Thomas interpolant on anisotropic elements. 
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1 Introduction 

In this paper, we study the Galerkin boundary element method (BEM) on graded meshes 
for numerical solution of the electric field integral equation (EFIE) on a Lipschitz polyhe¬ 
dral surface F in (i.e., F = dQ, where U C is a Lipschitz polyhedron). The EFIE 
models the scattering of time-harmonic electromagnetic waves at a perfect conductor, and 
the Galerkin BEM is widely used in engineering practice for simulation of this physical 
phenomenon. 

The Galerkin BEM considered in this paper employs divp-conforming lowest-order 
Raviart-Thomas surface elements to discretise the variational formulation of the EFIE 
(known as Rumsey’s principle). This approach is referred to as the natural BEM for the 
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EFIE (there exist other approaches, e.g., based on a stable mixed reformulation of Rum- 
sey’s principle, see [H]). Non-coercivity of the bilinear form in Rumsey’s principle (due 
to the inhnite-diniensional kernel of divr, cf. fl2.ip l signihcantly complicates the conver¬ 
gence analysis of Galerkin schemes. This problem can be overcome by using appropriate 
decompositions of vector helds in order to isolate the kernel of divp (we refer to discussion 
in [m Section 3], to an abstract theory in [T3l [12], and we outline available techniques 
for constructing such decompositions in Section H]). These ideas have led to major ad¬ 
vances in the convergence analysis and a priori error analysis of the BEM for the EFIE 
on (open and closed) Lipschitz surfaces, see fM[ ITBl [T^ [T^ [IS] for the h-version of the 
BEM and IHlIHllSllH] for high-order methods {p- and hp-BEM). All these results, however, 
assume shape-regularity of the underlying meshes on T. 

It is well-known that convergence rates of the h-BEM with quasi-uniform and shape- 
regular meshes are bounded by the poor regularity of solutions to the EFIE on non-smooth 
surfaces. For example, on a closed polyhedral surface T = dQ, the solution may be only 
H^(r)-regular (with a small e > 0 in the case of non-convex polyhedron G, cf. |20l Sec¬ 
tion 4.4.2]), and convergence rate of the h-BEM is only | -|- e in this case, whereas in the 
case of smooth solutions the lowest-order h-BEM converges with the optimal rate of | 
(see pn Theorem 8.2] and [S] Theorem 2.2]). Taking the cue from the h-BEM results for 
the Laplacian (see |28l [27]), we expect that an optimal convergence rate of the h-BEM for 
the EFIE can be recovered on the non-smooth surface T, if one employs the meshes that 
are appropriately graded towards the edges of T. These meshes contain highly anisotropic 
elements along the edges of T, and none of the results mentioned above is applicable in 
this case. Moreover, to the best of our knowledge, the quasi-optimality of the Galerkin 
h-BEM with graded meshes for the EFIE has not been studied in the literature, and with 
this paper we hll this theoretical gap. 

In the next section, we introduce necessary notation and formulate the EFIE in its 
variational form. In Section [3l we construct graded meshes on T, introduce the boundary 
element space, and formulate the main result of the paper— Theorem 13.11 — that estab¬ 
lishes quasi-optimal convergence of Galerkin solutions on graded meshes. The proof of 
Theorem 13.11 follows the approach suggested in |T2l US], summarised in [TSl Section 9.1], 
and extended to a general class of operators in [T^ Section 3]. At the heart of this ap¬ 
proach is the decomposition technique described in Section 0] Section [5] is instrumental in 
the construction of the corresponding discrete decomposition; here we establish new stabil¬ 
ity properties of the Raviart-Thomas interpolant of low-regular vector helds on anisotropic 
elements. In Section (6] we introduce the discrete decomposition and complete the proof 
of Theorem 13.11 An essential ingredient here is the projection operator Qh with enhanced 
approximation properties (see Proposition 16.Ij) . The proof of Proposition 16.11 is given in 
Section [71 
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2 The electric field integral equation 

The variational formulation of the EFIE is posed on the Hilbert space 
X = H-^/2^divr,r) := {u G divpu G 

Here, divp denotes the surface divergence operator, Hy is the dual space of H|'|'^^(r) 

(the tangential trace space of H^(r2) on T, see [TTl [T7]h and is the dual space of 

The space X is equipped with its graph norm || ■ ||x- We refer to o [la [m HE] 
for dehnitions and properties of H~^/^(divr, T) and other involved trace spaces. We also 
recall from min] that X is the natural tangential trace space of H(curl, H). 

In the present article, we use the same notation as in [5], where we recalled dehnitions of 
the full range of Sobolev spaces and differential operators needed for convergence analysis 
of the BEM for the EFIE (see Section 3.1 therein). In particular, we use a traditional 
notation for the Sobolev spaces (of scalar functions) (s G [—1,1]), Hq (s G (0,1]) 

and their norms on Lipschitz domains and surfaces (see [2^ ES])- The norm and inner 
product in LF‘{D) = H^{D) on a domain or surface D will be denoted by || • ||o,d and 
(•, •)o,_D, respectively. The notation (•, •)o,d will be used also for appropriate duality pairings 
extending the E^(El)-pairing for functions on D. 

For vector helds we will use boldface symbols (e.g., u = {ui,U 2 )), and the spaces (or 
sets) of vector helds are also denoted in boldface (e.g., H^(El) = with D C M^). 

The norms and inner products in these spaces are dehned componentwise. The notation for 
the Sobolev spaces of tangential vector helds on T follows OEsun]. In particular, L^(r) 
denotes the space of two-dimensional, tangential, square integrable vector helds on T. The 
norm and inner product in this space will be denoted by || • ||o,r and (•, •)o,r 5 respectively, 
and we will also use (•, Oo,? for appropriate duality pairings extending the (r)-pairing 
for tangential vector helds on T. The similarity of this notation with the one for scalar 
functions should not lead to any confusion, as the meaning will always be clear from the 
context . 

For a hxed wave number k > 0 and for a given source functional f G X', the variational 
formulation for the EFIE reads as: find a complex tangential field u G X such that 

a(u, v) := (Tfcdivr u, divr v) —/c^(^'fcU, v) = (f, v) Vv G X. (2.1) 

Here, (resp., T'fc) denote the scalar (resp., vectorial) single layer boundary integral op¬ 
erator on T for the Helmholtz operator — A —/c^, see [T6l Section 4.1] (resp., [T8l Section 5]). 

To ensure the uniqueness of the solution to fl2.ip we always assume that /c^ is not an 
electrical eigenvalue of the interior problem in H. 

3 Galerkin BEM on graded meshes. The main result. 

For approximate solution of (12. Ih we apply the natural BEM based on Galerkin discreti¬ 
sations with lowest-order Raviart-Thomas spaces on graded meshes. 
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First, let us describe the construction of graded meshes on individual faces of F. Here, 
we follow |281 Section 3]. For simplicity, we can assume that all faces of F are triangles. On 
general polygonal faces the construction is similar, or one can hrst subdivide the polygon 
into triangles. On a triangular face F C F, we hrst draw three lines through the centroid 
and parallel to the sides of F. This makes F divided into three parallelograms and three 
triangles (see Figured]). Each of the three parallelograms can be mapped onto the unit 
square Q = (0,1)^ by a linear transformation such that the vertex (0, 0) of Q is the image 
of a vertex of F. Analogously, each of the three sub-triangles can be mapped onto the unit 
triangle T = {x = {xi,X 2 )] 0 < Xi < 1, 0 < X 2 < Xi} C Q such that the vertex (1,1) 
of T is the image of the centroid of F. Next, the graded mesh on Q (and hence on T) is 
generated by the lines 

Here, > 1 is the grading parameter, and > 1 corresponds to the level of rehnement. 
Mapping each cell of these meshes back onto the face F, we obtain a graded mesh of 
triangles and parallelograms on F (see Figured])- Note that the diameter of the largest 
element of this mesh is proportional to Hence, h = 1/N dehnes the mesh parameter, 

and we will denote by T = {A^} a family of graded meshes = {iF; \JK = F} generated 
on F by following the procedure described above. 

Let us now introduce the boundary element space X/^. It is known that Raviart-Thomas 
surface elements provide an affine equivalent family of divr-conforming hnite elements 
under the Piola transformation, see [TTl Section HI.3]. We will write IVTq{K) for the local 
lowest-order Raviart-Thomas space on a generic (triangular or quadrilateral) element K, 
and we denote by = 7?.T"o(A^) the corresponding space of divr-conforming boundary 
elements over the graded mesh A^. 

The following theorem states the unique solvability and quasi-optimal convergence of 
the Galerkin BEM on graded meshes for the EFIE. 

Theorem 3.1 There exists ho < 1 such that for any f G X' and for any graded mesh A^ 
with h < ho and f3 G [1,3), the Galerkin boundary element discretisation of fl2.ip admits 
a unique solution Uh G X/, and the h-version of the Galerkin BEM on graded meshes A^ 
converges quasi-optimally, i.e., 

||u-u/,||x<G inf ||u-v||x, (3.1) 

vGX,i 

where the constant C may depend only on the geometry o/F and the grading parameter f3. 

The proof of Theorem 13.11 relies on an abstract theory for analysing convergence of 
Galerkin discretisations for non-coercive variational problems like (12.R . This theory was 
developed in [T3] . in, |T8l Section 9.1], and in [T^ Section 3]. In particular, it follows 
from the latter article that in order to prove Theorem 13. II we need to establish the following 
properties: 
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Figure 1: Graded mesh on the triangular face F C T. The triangular (resp., parallelogram) 
block of elements Tp (resp., Qp) is the image of the graded mesh on the unit triangle T 
(resp., the unit square Q). 

(A) the existence of a stable direct decomposition X = V 0 W such that a|vxv and 
— a|wxw are both X-coercive, and a|vxw and a|wxv are both compact; 

(B) the existence of the corresponding discrete decomposition X/^ = W/j, W/^ C W, 
that is uniformly stable with respect to the mesh parameter h; 

(C) the gap property 

llv — Vi, II Y 

sup inf -—77 -< s{h) with e{h) —)■ 0 as h —)■ 0. (3.2) 

VhGV,.veV ||V/i||x 

We will prove Theorem 13.II bv verifying these properties in Sections (Hand [6] below. 

Remark 3.1 Theorem 13.II remains valid ifT is a piecewise plane orientable open surface 
(see [13] for the problem formulation and the underlying tangential trace spaces in this case). 
The proof repeats the arguments in Sections^ and |6] below by using a specific construction 
of the decomposition X = V © W as described in [HI Section 3]. 

Remark 3.2 If some information about the regularity of the solution u to (12. Ij) is available, 
then convergence result of Theorem 13.11 translates into an a priori error estimate in the 
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natural 'K-norm. For scattering problems with sujfieiently smooth source functional f (e.g., 
with f representing the excitation by an incident plane wave), the regularity of the solution 
depends only on the geometry ofT. In particular, nonsmoothness ofT leads to singularities 
in the solution of the EFIE, severely affecting convergence rates of the h-BEM on shape- 
regular meshes. However, similar to the case of the Laplacian in [251127] . by employing 
the graded meshes with sufficiently large grading parameter (3 (depending on the strength of 
singularities in \i) one may hope to recover the optimal convergence rate (i.e., the rate of 
the h-BEM on guasi-uniform meshes in the case of a smooth solution). The main guestion 
here is whether the restriction on the grading parameter f3 that guarantees guasi-optimality 
of the Galerkin BEM in Theorem \‘d.l[ is sufficient for recovering this optimal convergence. 
We will address this issue in the forthcoming article jlO) . 

Throughout the paper, C, Ci, etc. denote generic positive constants that are indepen¬ 
dent of the mesh parameter h and involved functions but may depend on the geometry of 
r and the grading parameter (3. We will also write a < b and A B, which means the 
existence of generic positive constants C, Ci, C 2 such that a < Cb and CiB < A < C 2 B, 
respectively. 


4 Decomposition technique 

Let us address property (A) in Section [3l One way to obtain a suitable decomposition is 
to employ the (r)-orthogonal Hodge decomposition of X, cf. [I5]. In the context of the 
h-BEM on shape-regular meshes, this idea was successfully exploited in [T6l 1211 [191 113] • 
However, in the case of non-smooth snrfaces, the regnlarity of surface gradients in the 
V-component of the decomposition may be poor, and this causes substantial technical 
difficulties in the analysis of the p- and hp-BEM. For the p-BEM on plane open screens, 
a modihcation of the above strategy was suggested in |6], where we consistently used 
the -inner product and proved H ^/^-orthogonality of the Hodge decomposition. 

Unfortunately, these ideas do not generalise immediately to closed polyhedral surfaces or 
piecewise plane open screens, neither to the hp-BEM with quasi-uniform meshes, cf. j^. 
When attempting to use the Hodge decomposition for convergence analysis of the h-BEM 
on graded meshes, poor regularity of vector helds in the V-component leads to severe 
restrictions on the grading parameter (3. 

An alternative technique employs a regularising projection R : X —)■ X to construct a 
decomposition of X with enhanced regnlarity of the V-component (see |23| . [T8l Section 3], 
and m Section 4.3.1]). The projection R is dehned by employing the H^(r2)-regular vector 
potentials from the following lemma. 

Lemma 4.1 [21 Section 3] For any bounded Lipschitz domain hi C there exists a 
continuous mapping L : curlH(curl, hi) —)■ H^(r2) such that curl L $ = $ for all $ G 

curlH(curl, H). 
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Let u G X. The regularising projection R : X i—)■ X is defined as follows: 

Ru := ((L$) X i^)|r, 

where L is from Lemma 14.11 denotes the unit outward normal to hi, $ := Vtc, and 
w G H^{Q) is the solution to the problem 


—Aw = 0 in n, 

Vw ■ u = divru on L. 

The fact that divpu dS = 0 for each connected component S of T guarantees $ G 
curl H(curl, n), and Lemma [4.11 can be applied. By using elliptic lifting theorems, trace 
theorems, and the continuity of L, we conclude: 

3C' = C'(r)>0 such that || Ru||jji/ 2 ^p^ < C ||divru||j:^-i/ 2 (r) VugX, (4.1) 

where Hy^(r) C X is the rotated tangential trace space of H^(r2) on T = dfl, see [TTl [T7] . 
By construction of R, we have on T 

divr Ru = divru for any u G X, (4-2) 


that is R^ = R. 

Now we can define the decomposition 

X = V © W with V := R(X) C Hy^(r) and W := (Id - R)(X). (4.3) 

Decomposition fl4.3p was used in |9] to prove unique solvability and quasi-optimal conver¬ 
gence of the hp-BEM with locally variable polynomial degrees on shape-regular meshes. 
As we will see in this paper, the same decomposition technique can be used effectively in 
the analysis of the h-BEM on graded meshes. 

By fl4.2p we conclude that W comprises divp-free vector fields. Stability of decompo- 
sition (14.3p follows from inequality (14.ip and the continuous embedding Hp (T) ^ X. 
Furthermore, the embedding V L^(r) is compact by (14. ip and Rellich’s theorem. Thus, 
thanks to the iL“^/^(r)-coercivity (resp., H_|_^'^^(r)-coercivity) of 4/^ (resp. ^k) (see [T8l 
Lemmas 8, 7]), the X-coercivity of a|vxv and a|wxw is proved by the same arguments as 
in [T^ proof of Theorem 3.4]. The compactness of a|vxw and a|wxv is due to the conti¬ 
nuity of : Hj)^'^^(r) —)■ Hy^(r) and the compactness of the embedding V ^ Hj)^'^^(r) 
(see [ISl Lemma 9]). This proves (A). 

Before we can define the discrete counterpart of decomposition (14.3p . we need to find 
a suitable projector onto the space of Raviart-Thomas surface elements. Stability of the 
discrete decomposition will follow from stability properties of the Raviart-Thomas inter¬ 
polation on anisotropic elements, which is the subject of the next section. 


7 













5 Raviart-Thomas interpolation on anisotropic 
elements 


In this section, we establish new stability properties of the Raviart-Thonias interpolant on 
anisotropic elements. We will also prove the corresponding interpolation error estimates. 

In the context of the hnite element method, the analysis of interpolation operators on 
anisotropic elements can be found in EHH. For the Raviart-Thomas interpolation (see, 
e.g., [H Section 3]), the main idea is to study componentwise stability of the interpolant 
on a reference element K. For sufficiently regular vector helds, this study relies on the 
fact that the standard (scalar) trace operator is well dehned for functions in IF^’^(iF) for 
any p > 1. In our BEM application, however, the stability result is needed for low-regular 
vector helds living in fractional Sobolev spaces H^(iF) ft H(div,iF) with 0 < s < 1/2. 
For such vector helds, the trace of the normal component only exists in a weak sense. 
Therefore, instead of the standard trace argument used in |T], we use Green’s formula (15.3p 
and consistently employ the anisotropic seminorms dehned below. More precisely, for 
s G (0,1/2] we introduce the Ff^-seminorms of anisotropic type. On the reference square 
Q = (0,1)^ these are dehned as follows: 


~ OIr/yo,!)'^^1- 

0 


These dehnitions are meaningful for all u G H^{Q) due to [25| Theorem 10.2] which yields 
that 

On the reference triangle T = {{xi, X 2 ); 0 < xi < 1, 0 < X 2 < xi}, the following seminorms 


u 


2 

AHlif) 


m 


•,^2)1 




U 


2 

AH0) 


■) lHqO,a:i) dXi 


are also well dehned for all u G H^{T). Indeed, by Theorem 1.4.3.1 
continuous linear operator (called the extension operator) E ; H^{T) 


of [21] there exists a 
—)■ H^(Q) such that 


Eu\f = u, WueH^f). 


Hence, applying (15.ip to Em, we have 





We recall that the Raviart-Thomas interpolant IIrtu is well-defined for any u G H^(R') 
(s > 0) snch that divu G L‘^{K), where K is any triangle or rectangle. Indeed, for snch 
vector fields the following Green’s formula has a meaning (see, e.g., [3 Lemma 2 . 1 ]) 

{u,\/(p)o,K+ divu(^ = (u ■ n, v9)o,aK, e (5.3) 

Jk 

with e G (0, s) and n denoting the outward unit normal to dK. Hence, taking ip G 

such that 99 = 1 on the edge e C dK and 99 = 0 on dK\e, we can define (u ■ n, l)o,e := 

(u ■ n, (^)o,ax- 

In what follows, we will denote by Hrt the Raviart-Thomas interpolation operator on 
the reference element K {K = Q or T). 

5.1 The reference square 

On the reference square Q, we denote by fii and 63 the edges parallel to the a:i-axis, and 
by 62 and 64 the edges parallel to the a: 2 -axis (see Figure [2]). We recall from [11] that the 
lowest order Raviart-Thomas elements on Q are defined as 

o{Q) = {{a + cxi,h + dx 2 Y] a,5,c,d G M} , 

and that the associated degrees of freedom are given by 



Figure 2; The reference square Q and the reference triangle T. 
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Theorem 5.1 (i) For all u G H®(Q) with s > 1/2, we have for / = 1, 2 


||(nRTu)z||oQ < (5.4) 

(ii) For all u G H'^((5) H H(div,(5) with 0 < s < 1/2, we have for I = 1, 2 

||(nRTu)z||oQ < WuiWhs^q) + + ||divu||o^^ (mod. 2). (5.5) 

Proof. By symmetry, it will suffice to prove both statements for / = 1. One has 

IIrtu = (a + cxi, b + dx 2 )~^ with a, b,c,dG M. 


Hence 

||(nRTu)i||oQ = \\a + cxi\\QQ < \a\ + |c|, 
and it remains to estimate |a| and |c|. Observe that 

a = -(u • n, l)o,e^, and a + c = (u • n, l)o,g^ 2 . 

Again by symmetry, it will suffice to estimate (u ■ n, 1 ) 0 , 64 . 

If u G H'^(Q) with s > 1/2, then the trace of Ui on 64 is well-defined, and 

l®l ~ |(u ■ n, 1)0,64! ~ 1)0,64! ~ I!'*^il!ir®-i/ 2 (e 4 ) ~ l!’*^il!iy'>(Q)- 

This proves statement (i) (for Z = 1). 

Now, let us consider u G H'^(Q) O H(div,Q) for 0 < s < 1/2. In order to estimate 
(u ■ n, 1 ) 0,64 ffi this case, we fix a function (p G i/^“^(Q), e G (0, s), such that 99 = 1 on 64 
and 99 = 0 on dQ\ei 4 . Then by Green’s formula fl5.3|) we have 


(u-n, 1 ) 0,64 = (u,V(^)oq-F ldivuip = {uudiip)f^Q + {u 2 ,d 2 ip)QQ+ divuip. (5.6) 


'Q 


'Q 


For the term {ui,diip)QQ, we first use a standard duality argument, 

!(wi, (9l<y9)Q q! < |!ni||^e(Q) |!<9l<y9|!j^_e(Q), 

and by the continuity property of di : H~^{Q) (see Theorem 1.4.4.6 of |21| ) 

we find 

!(ui, 51^9)0 q! < I!'*^iI!_h's(q) y^di!( 5 - 7 ) 

The term (M2,f^2<^)oQ requires more subtle analysis. First by (15.11) and Fubini’s theorem, 
we can write (hereafter, I = (0,1)) 


{U2, - 92 <^)o,q = 


(^ 2 ( 3 : 1 , ■),d 2 (p{xi, ■))Qjdxi. 


(5.8) 
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Now we use a density argument to show that 


d2(p{xi, ■))Q fdxi = 0, 


(5.9) 


for all h G Indeed, £x a sequence of smooth function ^pn such that 

in as n —)> oo. 

Then for all Xi E I = (0,1), we have 

{l,d2(Pn{Xl,-))oT= / d2(Pn{Xl,X2)dX2 = (Pn{Xl,l)-^n{Xl,0). 

Jo 

For any h E multiplying this identity by h and integrating the result in Xi E (0,1), 

we obtain ^ 

/ h{Xi){l,d2(Pn{Xl,-))ojdXi = / h(pn- / h(pn. 

Jo J J 

Hence, as n —)■ oo we find that 


h{xi){l,d2(p{xi,-))Qfdxi = I hip- I hip, 


which proves fl5.9p by recalling that 99 = 0 on hi and 63 . 
Coming back to fl5.8p and using fl5.9p . we have 


{u2,d2ip)QQ= / {U2{xi,-) - J^j{u2{Xi,-)),d2ip{Xi,-))^Y^^i^ 

Jo 

where JAf{u 2 {xi, ■)) = U 2 (xi, X 2 ) dx 2 is the mean of U 2 (xi,-) on I = (0,1) (clearly 
J'^f(u2(xi, ■)) E At this stage we hrst use a duality argument and then again 

Theorem 1.4.4 .6 of [ 2 T] to obtain 


i(u2,d2(p) 


0,qI - 


\\u 2 {xi, ■) - J^y{u 2 {xi, •)||nqo,i)||52<^(a^i, Olln-qo,!) dxi 


< / - ■MT{u2{^iE)\\H^{o,i)\\^{xi,-)\\H^-^{o,i)dxi. (5.10) 

Jo 

We use Friedrichs’ inequality to estimate 

\\u 2 {xi, ■) - J^j{u 2 {xi, Olinqo,!) < 11^2(2:1, •) - J^y{u 2 {xi, Olinqo,!) < \u 2 {xi, ■)\h‘>{o,i)- 
Using this estimate in (I5.10p and applying the Cauchy-Schwarz inequality we arrive at 


|(m2,'92<^)o,q| < |M2Ujy|(Q)( / ||<^(xi, •) dx 


(5.11) 
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The last term on the right-hand side of fl5.6p is estimated by applying the Cauchy-Schwarz 
inequality: 

divuLp < ||divu||o_^||(p||o^^. 

J Q 

Using this estimate and inequalities (15.7p . (Ih.llh in (Ih.hh . we obtain (I5.5h (for / = 1). □ 

Corollary 5.1 (i) For all u G H^(Q) with s > 1/2, we have for Z = 1, 2 

\Wl - (nRTu)zllog < \ui\jjs(^Qy 

(ii) For all u G H®((5) H H(div, Q) with 0 < s < 1/2, we have for / = 1, 2 

\\ui - (nRTu);||o^^ < + |Mi+i|^^.^^(Q) + ||divu||o^^ (mod. 2). 

Proof. It is sufficient to prove only statement (ii). For / = 1, we take u = u — (A1 
where M.qUi = JqUi. One has 

u — IIrtu = u — IIrtu, div u = div u, 

and 

I“iIh1 /2(Q) = |hl|^l/2 (Q) ^ ||'Ui||j^1 /2 (q). 

The assertion then follows by applying estimate fl5.5p to u. The proof is analogous for 
1 = 2 . □ 

Corollary 5.2 Let s G (0,1/2]. For all u G H^(Q) such that divu G M, we have for 

1 = 1,2 

II(nRTu)illQQ ^ -f (mod. 2) (5-12) 

and 

\\ui — (nRTu);||g^^ < \'^i\h‘{Q) (j^od. 2). (5.13) 

Proof. Note that the function (p G Fl^~^{Q) in the proof of Theorem 15.11 (ii) can be 
chosen to have zero average on Q, i.e., /g<p = 0. Indeed, if this is not the case then we 
£x G C^{Q) such that = 1, and consider (p — (/^(p)'0 G that has the 

same values on dQ as ip and also has zero average on Q. Then, for divu G M, the last 
term on the right-hand side of fl5.6p vanishes, and fl5.12p (resp., fl5.13p l follows by the same 
arguments as in the proof of Theorem 15.11 (ii) (resp.. Corollary 15.11 (h)). □ 
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5.2 The reference triangle 

On the references triangle T, we denote by ei (resp. 62 ) the edge on the xi-axis (resp., 
parallel to the X 2 -axis), and by 63 the oblique edge (see Figure [ 2 ]). We recall from [TT] that 
the lowest order Raviart-Thomas elements on T are dehned as 

IZToiT) = {{a,by + c(a;i,a; 2 )'^; a, 6 ,c G M} , 

and that the associated degrees of freedom are given by 

/ u ■ nds, z = 1, 2, 3. 

Jci 

Theorem 5.2 (i) For all u G H'^(T) fl H(div, T) with s > 1/2, we have for / = 1, 2 

||(nRTu)/||of < + ||divu||o_^. (5.14) 

(ii) For all u G H'^(T) 0 H(div,T) with 0 < s < 1/2, we have for / = 1, 2 

||(nRTu)i||of < + l|divu||o^^ (mod. 2). (5.15) 

Proof. We will prove both statements for / = 1. The proof is analogous in the case 1 = 2. 
One has IIrtu = (o, b)~^ + c{xi, X 2 Y with a,b,c ^ M. Hence 

||(nRTu)i||o^f = \\a + cxillo^^ < |a| + |c| < |a + c| + 2 |c|. 


Observe that 


2c = div Hrtu 


By the Cauchy-Schwarz inequality this yields 


2 


J div u. 


|c| < ||divu||o_^. 

Therefore in order to estimate ||(nRTu)i||g^ one needs to bound |a + c| = |(u ■ n, l)o,e 2 l- 
If u G H®(T) with s > 1/2, then we can use trace theorem to estimate 

|(u-n, 1)0^62 1 I(^1) 1)0,62 1 ^ ll"^!IIirs(T)■ 

This proves statement (i) (for / = 1). 

Now, let us consider u G H'^(T) O H(div, T) for 0 < s < 1/2. In order to estimate 
(u ■ n, 1 ) 0,62 in this case, we £x a function (p G F[^~^{T), e G (0, s), such that (p = 1 on 62 
and (p = 0 on dT\e 2 - Then by Green’s formula (15.31) we have 

(u • n, l)o ,^2 = (m, V<^)of + / divu(y9 = (Mi,ai<p)o^ + (m2,52</2)o^ + / divu(p. (5.16) 

Jf ’ ’ Jf 
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As before (cf. fl5.7p ). it is easy to show that 


The term (M2,c?2V^)of treated as in the case of the reference sqnare except that the 
interval I = (0,1) in the a; 2 -variable will be replaced by the interval I{xi) := (0,xi), but 
the arguments remain mainly unchanged because xi < 1. More precisely, hrst by (15.21) and 
Fubini’s theorem, we may write 


{u2,d2^p)Qf= / {u2{xi,-),d2(p{xi,-))oj(xi)dxi. 


(5.18) 


With a property similar to (15.9h we deduce that 

{u 2 ,d 2 Cp)o^f = {U 2 {xi,-) - Mi(xi){u 2 {xi, ■)), d 2 cp{xi, dxi, 

where Aii(^xi){u 2 {xi, ■)) = ^ U 2 (xi, X 2 ) dx 2 is the mean of U 2 (xi, ■) on I(xi) = (0,xi). 
Then a duality argument yields 

j(u2,d2cp)gfj< / ||M2(a:i, •) - Al/(xi)(M2(a;i, •))lkdh^i))ll^2</?(^i’Olln-q/Cxi)) dxi. (5.19) 
Jo 

Using a scaling argument and Friedrichs’ inequality we estimate 


\\U2{XI, •) - J^I{xi){u2{Xi, •))||Hq/(xi)) < \\U2{XI, •) - J^I[xi){u2{Xi, ') ) II Hq/(xi)) 

< C\u2{xi, (5.20) 


with C > 0 independent of xi G (0,1). 

For the second factor in the integrand in fl5.19p . in order to apply Theorem 1.4.4.6 of 
[21] on a hxed domain, we hrst notice that Lp{xi,-) E IIq~^(I(xi)) a. e. in (0,1) 9 Xi. 
Therefore, for almost all xi E (0,1) there exists a sequence of functions cpn ^ C'“(/(a:i)), 
n = 1,2 ,..., such that cpn ^ 9 ^(xi, ■) in FrQ“^(/(a;i)) as n —)■ 00 and 

d(Pn(x 2 ) dx 2 = 0 Vn = l,2,.... 

Using a scaling argument, we deduce that 



\\d^n\\ 


d^Pn V dX2 

= sup - 

vp^O 


^ ~2 

< Xi ^ sup 

vj^O 


/q dipn V dX 2 

I'^|h®(0,1) 


As d(pn has zero average, we can estimate 


||^V^n||_ff-®(/(a:i)) < X^ 




d^nVdX2 


sup 


^^0- fo ^^=0 


< 7 . 2 += 


/q d(pnVdX2 


sup 

v&H^(o,i) ITI|irs(o,i) 


<X]^ 


\\d(p 


n II • 


^^0, /o v=0 
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Hence, by Theorem 1.4.4.6 of |21) we prove that 

--+e ^ --+e ^ 

||^V^n||i/-^(7(xi)) ^ ^ ||</2n||//l-e(0,l) < ^ | (p„ | j^l-e( 04 ). 

Mapping back to the interval /(xi) = (0,xi) we have 


As n —)■ cxD we hnd 

\\d 2 ip{xi, < Ci|(^(xi, ■)\m-'^{i{x^)) a. e. on (0,1) 9 xi, 

with Cl > 0 independent of Xi. 

Using estimates fl5.2Up and fl5.2ip in fl5.19p we arrive at 


( 6 , 21 ) 


\(U2,d2lfi) 


0,T| ~ 


|W2(X1, •)|fl'hO,xi)||‘^(xi, OIIhi-TD.^h) '^^1- 


Then the Canchy-Schwarz inequality yields 


\{U2,d2ip)Q^f\ ~ l“2Ui^|(f)( / y{xi,-)\\m-e(o,x,)dx 


(5.22) 


The last term on the right-hand side of fl5.16p is estimated by using the Cauchy-Schwarz 
inequality: 

divuif < ||divu||o^^||9?||o56. 


Using this estimate and inequalities (I5.17p . (I5.22p in fl5.16p . we obtain fl5.15p (for /=1). □ 


Corollary 5.3 (i) For all u G H'^(T) fl H(div, T) with s > 1/2, we have for I = 1, 2 

\\ui - (nRTu);||o_^ < + ||divu||o^^. 

(ii) For all u G H^(T) fl H(div, T) with 0 < s < 1/2, we have for Z = 1, 2 

\\ui - (nRTu);||o^^ < \ui\j^s^f^ + + ||divu||o56 (mod. 2 ). 

The proof of this statement is similar to the proof of Corollary 15.11 

Counterexample 5.1 Here we provide a counterexample which demonstrates that for low- 
regular vector fields the terms \ui+i\^jjs (q) in fl5.5p and \ui+i\^jjs (f) in fl5.15p cannot be 
omitted. This is in contrast to the case of sufficiently-regular fields in fl5.4D and fl5.14p (see 
also Lemma 3.3 in [1]/. In particular, if we assume that 

||(nRTu)2||of < 11 M2 111^1/2 (f) + ||divu||o^^ (5.23) 
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for all u G H 2 (T) fl H(div,T), then we will arrive at a contradiction. Indeed, inspired by 
Example 2.6 of [3], we define onT x (0,1) 

v^{xi,X 2 ,X 3 ) = {xi — l)w‘^{x 2 ,X 3 ) with Xs) = min 11, £ loglog -1 for any e > 0. 

Here, r = (x^+x"^)^, and e is the Euler number. Taking Wv’^xv onT (here v = (0, 0, — l)j, 
we find a divergence-free vector field 

u^(xi,a:2) = ((1 -Xi)a2W^(x2,0),M;^(x2,0))"^. 

Simple calculations show that = (0,1)"*", and by trace theorem we have 

Since x(o i)) —^ 0 as e —)■ 0 (see j3l Example 2.6]j, we conclude that (u ^)2 —^ 0 

in H^lT) as e ^ 0. This seems to contradict (I5.23p but not directly because the first 
component of is not in H 2 (T). Hence, in order to arrive at a contradiction, we need 
to show that if fl5.23p holds for all u G (T) n H(div, T), then satisfies (I5.23P (with a 
constant indepedent of e). Indeed, for a fixed e > 0, as 0) G H^{f), 1), we can consider 
a seguence of smooth functions Wn G C°°([0, 1]) such that 

Wn w^{-, 0) in H^{0, 1) as n ^ 00 . 


Then we define 

Un{Xi,X2) = {{I - Xi)dWn{x2),Wnix2)y. 

One has u„ G H 2 (T) and divu„ = 0. Moreover, IlRTUn — )■ as n ^ 00 . Therefore, 

applying estimate fl5.23p to and letting n —)■ 00 , we conclude that satisfies fl5.23p . 

Remark 5.1 By Counterexample^^ we can easily show that a result similar to Lemma 3.3 
in jl] for the Nedelec interpolant on the tetrahedron T 3 = {(xi, 0 : 2 , 0 : 3 ) G Xj > 0, i = 
1, 2, 3 and 0 < Xi + X 2 + X 3 < 1} is not valid. In other words, the anisotropic estimate 

II (nNedV)zllQ ^ ll'^dlHl(T3) II II 0,f3 ’ / = 1, 2, 3 

does not hold for all v G H^(T 3 ) = (^H^{T^)Y■ 

5.3 Anisotropic elements 

In this subsection we will denote the functions on the elements K and 7^ by u and u, 
respectively. Analogous notation will be used for coordinates (e.g., x G TT and x G K) and 
for differential operators (e.g., div and div). 

First, let us prove the following auxiliary result. 
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Lemma 5.1 Let K be the image of the reference element K (K = T or K = Q) under 
diagonal scaling with matrix B = , where hi > 0. Then for any u G there 

holds 

max{hY\hl^^} 




hih2 


where u(x) = u(i?x), x = {xi,X 2 ) G K. 
Proof. By the definition of the dnal norm 


\u\ 


H-i/ 2 (^) - sup 




ll^ll 

v£H^/^{K) ll'^ll_ffl/2(X) 


(5.25) 


We now estimate the norm ||n||^i/ 2 (-|f). If n G H^{K), then diagonal scaling yields 


m = ihih2)-^ \\v"^ 


' 0 ,K 

,-i no „,ii2 


o,Ky 


ll-9iu|lo,x - ^1^2 II-9 iu||o,k, l|-92u|lo,i? - ^1 ^2 11-92^115,;^. 


Therefore, 






' 0 ,K 


and by interpolation between and we find that 


1-112 min{hi,h2} 2 

|U||^V2(^)^ --11^1' 


fi\h2 


|U||2 


HhlKy 


'iv e (5.26) 


h^h2 max{hi,h2}" 

Since {u,v)q^ = {hih 2 )~^ {u,v)o^k, we nse fl5.26p in fl5.25p to obtain inequality fl5.24p . □ 

Now, we are in a position to prove the stability result and the corresponding error 
estimate for the Raviart-Thomas interpolation on anisotropic elements. 


Theorem 5.3 Let K he either the triangle T with vertices (0,0), (hi,0), {hi,h 2 ), or the 
rectangle Q with vertices (0,0), (hi,0), (0, h. 2 ), {hi,h 2 ), where hi > 0. Denote h^ax ■ = 
max {hi, h 2 }- Then for any u G H^/^(R') with divu G M there holds for I = 1, 2 


||(I1rtu);||q i^ < 


0 ,K 


+ 


h^ 

max 

hih2 


+ ||divu||^. 


l/2(i^) 


and 


(mod. 2) 
(5.27) 


\ui — (IIrtu) 


|2 

lo,;^ 


< 


h^ 

max 

hi/i2 


|Ui|^i/2(;^) + |u/+l| 


AHllliK) 


+ l|divu||^_i/ 2 (^)^ (mod. 2). 

(5.28) 
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Proof. We consider only the case of the triangle, K = T. The proof is similar in the case 
K = Q. 

We nse the Piola transformation to dehne u E fl H(div, T) on the reference 

triangle T as follows: 


u(x) = hih 2 B ^u(i?x) with B = 


hi 0 

0 h2 


Then we have 


\Ui\ 


0,T 


— /i 2 (hi/12) ||'*^l|lo,T — ^2 ll'^i|lo,r) 





IT JT 

,2 tu . ^-2 



|mi(: 


wi(y)P 


(1^1 - 

-yi\ 

|2 + 

\X2 - y2\ 

2)3/2 


= h2{hih2y 



dxdy 

|mi(x) -ni(y )|2 


tJt (hi Vi - 2/iP + ^2 V 2 - 


-dxdy 


and 


|2 _ 

r,: 


Jo ' 


— hi \Ui\jjl/2(j^y 




= hi 


-1 


rhi 


r-haxi/hi i‘h2Xi/hi /i2|^2(Xi,X 2) - M2(a;i, 2 / 2 ) P dX2 dy2 


h2 V 2 - 2 / 2 I 


hi 


dxi 


hi \u2y 


AH^yT) 


Fnrthermore, the standard properties of the Piola transformation yield 


nRTu(x) = ——5nRTu(x) = 

hih2 


l/ha 0 

0 1/hi 


nRTu(x) and div u = hih 2 div u G 


Therefore, applying Theorem 15.21 (ii) and Lemma 15.11 and nsing the fact that div u G M 
(hence, ||divu|V — ll^ivu||i:i_i/ 2 (f)), we obtain 


||(nRTu)i||o,T “ h2 hih2 ||(nFi,Tu)i||Q 


< hih 


-1 






+ + lldivu 


< 


|2 


if 


I ||2 I h^a^ I |2 I hi I |2 I hi , II ,. ||2 


Recalling that hmax = max {hi, h 2 }, it is easy to see that 


ho ho 


< 


h^ 

' ‘'max 
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and then ineqnality (15.271) follows (for I = 1 and K = T). 

Arguing as above and using Corollary 15.31 (ii) instead of Theorem 15.21 (ii) we establish 
the error estimate in (15.281) for / = 1. The proof is analogous in the case 1 = 2. □ 


We can now estimate the L^-error of the Raviart-Thomas interpolation on the graded 
mesh on T. The specihc estimate that we need is for H^/^-regular vector helds with 
discrete divergence. 

1 /2 

Lemma 5.2 For any u G (T) such that divpu G divrX/i there holds 

||u - IlRTullo.r < ||u||jjy 2 (p). (5.29) 

Proof. Let F be a face of T, and let Tp G F he a triangular block of elements, see Figured] 
(the arguments are analogous for the parallelogram block of elements Qp). The triangle 
Tp is mapped onto the unit triangle T by the affine transformation which is independent 
of the mesh parameter h. Let us hrst establish the error estimate for the Raviart-Thomas 
interpolation on the unit triangle T partitioned into elements as shown in Figured! 

The graded mesh on T comprises the quadrilaterals Kij = = i > j) 

isomorphic to (0, hi) x (0, hj) with h^ > hj and the triangles Ka isomorphic to the triangle 
with vertices (0,0), (0,hj), {hi, hi). Applying error estimates from Theorem 15.31 on each 
element Kij {i > j), we have for / = 1, 2: 

\\ui - < h^,hj^ + ||div(mod. 2). 

Summing these estimates over all elements in T and recalling that hfhj^ < for 

1 A LJ A we obtain 


u 


hlRTull 


2 ^ 

0,T ~ 


TV . 2 

i,j=l ^ l=l 

i>j 




+ ||divu||^_i/ 2 (^^.) ). (5.30) 


Note that by standard superposition argument 


N i ^ r 

^ E/i ■^2(3^1, ■) Ih 1 / 2 (o,xi) 

i=l j=l i=l 

/•I 


ll531l 


'0 


\U2{Xl,-)\jjl/2^Q^^^^dXl — \U2\^^l/2^.f^ < ll'“2||j:^l/2(f), (5.31) 


and similarly for ui. Hence, using standard superadditivity properties of the H^^^-semi- 
norm and the iL“^/^-norm, we deduce from (I5.30p the following error estimate on T; 


|u - HrtuIIo^^ < h^ (^||u||jji/ 2 (f) + ||divu||^_i/ 2 (^)) 
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Applying now the Piola transformation associated with the mapping Tp —)■ T, patching 
together all individual blocks of elements on all faces of F, and using the superadditivity 
of if ^- 1/2 -norms (as the functions of subdomains), we obtain 

||u - nRTu||o,r < (l|u|lHi/"(r) + l|divru||p-i/ 2 (p) j (5.32) 

(here and below we use the space Hl.(r) , s > 0, which is dehned in a piecewise fashion by 
localisation to each face of F, with the norm ||u||^s := ll^li^llH®(F))- 

Inequality (I5.29p follows from fl5.32p due to the continuity property of divp: Hy^(r) 
^-i/2(r) ("ggg Section 4.2]). □ 


6 Discrete decomposition and the gap property 

Following the ideas from [13] and [T8l Section 9.1], we can use the Raviart-Thonias inter¬ 
polation operator IIrt to dehne the discrete counterparts of V and W in fl4.3p (e.g., we 
can set V^ '■= nRT(R(X/i)), where R is the regularised projector introduced in Section 0]). 
However, as it follows from the results in Section [S] the Raviart-Thomas interpolation of 
low-regular vector helds on graded meshes is only stable (with respect to the L^-norm) 
when < 2. Since the dehnition of the energy space X for the EFIE involves the dual 

_1 /o 

space Hj| ' (F) with a weaker norm than ]| • ]|o,r, we can relax the restriction on the grad¬ 
ing parameter f3 by employing a different projection onto the boundary element space and 
using a duality argument on individual faces of F. This approach was successfully used by 
Hiptmair and Schwab in [241 Section 8] and by Buffa and Christiansen in [I3l Section 4.2.2] 
in the context of the h-BEM with shape-regular meshes for the EFIE (see mu for appli¬ 
cations of these ideas to the analysis of the p-BEM and the hp-BEM with quasi-uniform 
meshes). We will demonstrate below that using these techniques together with stability 
properties and error estimates for the Raviart-Thomas interpolation on anisotropic ele¬ 
ments, one can design a stable discrete decomposition of the boundary element space on 
and prove the corresponding gap property fl3.2p for any /3 < 3. 

The construction of the desired projection operator is technically involved. Therefore, 
we formulate here the hnal result relevant to our discussion and give a detailed proof in 
the next section. In the Froposition below, Hq denotes the L^(F)-projection onto the space 
of piecewise constant functions over the mesh A^, and H_^'^^(F) denotes the dual space of 
hF2(f) L^(F) as pivot space). 

Proposition 6.1 There exists an operator Qh : Hi(F) nH(divr,F) (s > 0) such 

that 

divr o Q/j = Ho o divp, (6.1) 

and for any e > 0 

||u - - nRTu]|H(divr,r) Vu G H1(F) fl H(divr, F). (6.2) 
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Thus, the operator Qh inherits the crucial commuting diagram property fl6.1l) of the 
classical RT-interpolation operator and, at the same time, allows to gain an extra power 
of h when estimating the error (u — Q/^u) in the dual norm. 

1 /o 

Corollary 6. 1 For any u G H_|_' (T) such that divru G divrX/j, one has G X/^, 
divrQ/iU = divrU, and for any £ > 0 there holds 

||u - Qhu||x = ||u - (6.3) 

Since || ■ |lH-i/2(r) ~ II ' llH"^/^(r)’ result immediately follows from Proposition 16.11 
and Lemma [5.21 due to the commuting diagram property for IIrt- 

Since RX/i C Hy^(r) (see (14.1|) ) and divpRX/j = divr X/^ (see (14.2|) ). the following 
definitions are valid thanks to Proposition 16.11 

V;, := {Qh o R)X;„ W;, := (Id - Qh o R)X^. 

Using the commuting diagram property fl6.ip . we have 

divrQh Ru.h = no(divr R uh) ^ no(divrU/j) = divrU/^ Vu/^ G X^. (6.4) 


Therefore, 


R Q/i R = R on X/i, 


(6.5) 


and hence o R : X/j —)■ X/^ is a projection. This fact conhrms that X/j = Vh © W/j. 
Property (16.41) also implies Wh C W, and Corollary 16.11 yields stability of the discrete 
decomposition in the following sense: there exists C = C{r,(3) such that 


lEDi 


Q/i R u/i||x < C* II R u/i|| jji/2,p. < C ||divrU/i||p^-i/2(-r) < C ||u/i||x Vu/^ G X/^, 


provided that (d <3. This verifies property (B) from Section [3l 

It remains to establish the gap property (C). Inequality (13.2p in (C) is another conse¬ 
quence of Corollary 16.11 

• r 11''" “ '^hWx , II R V/i — V/i||x IMl II Rv/i — Qh R V/i||x 

sup mf -77- < sup -- = SUp -- 

VhGV^vev ||v/,||x v^eVh l|vft||x v^ev^ \\^h\\x 

R V/i||„l/2.„. I |4.1t 

©) < f^3l2-yi2-e_ 


< sup 

v^GV^ 


V/i X 


This completes the proof of Theorem 13.11 
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7 Proof of Proposition 16.1 


In this section, we give a constructive proof of Proposition 16.11 For any u G Hi(r) fl 
H(divr, F) we construct QhVl in the Raviart-Thonias spaces on individual faces of F. Let F 
be a single face of F, and let A^(F) denote the restriction of the graded mesh onto F. 
For the sake of simplicity of notation we will omit the subscript F for differential operators 
over this face, e.g., we will write div for div^T’. We will also write (•, •) for the LP‘{F)- and 
L^(F)-inner products, and similarly || • || for the corresponding norms of scalar functions 
and vector helds. First, let us prove the following auxiliary result. 


Lemma 7.1 For any s > 1/2, the Raviart-Thomas interpolation operator IIrt : H^(F) fl 
H(div,F) —)■ TZT~q{A^{F)), isI-P{F)-stable, i.e., there exists a constant C > 0 independent 
of h such that 

||nRTu||o,i. < C'(||u||h»(^) + lldiv u||o,f) Vu G W{F) n H(div, F). (7.1) 


Proof. Similarly to the proof of Lemma 15.21 it is sufficient to establish fl7.ip for the unit 
triangle T partitioned into elements as shown in Figured] (this is because the affine trans¬ 
formations that map triangular blocks of elements Tp <zF onto T are independent of h). 

The graded mesh on T (see Figured]) comprises anisotropic quadrilaterals Kij = /* x Ij 
{i,j = i > j) isomorphic to (0, h*) x {0,hj) with hi > hj and shape-regular 

triangles Ka isomorphic to the triangle with vertices (0,0), (0,hj), [hi^hf). Using the 
Piola transform associated with the mapping K —)■ K^j {i > j; K = T or Q), we dehne 
u G W{K) n H(div, K). Then, by the standard properties of the Piola transform we have 

ll^illoA - WMlxip l|divu||2^^ ~ hihj lldivu||2_^... 


||(nRTu)i||o,^.. ~ hi hj ^ IKUrtG) 
The application of the scaling argument yields; 


'0,K' 


Ml 


2 

Lf/qf) 




l“l( 

X)- 

w(y)P 


(1^1 - 

-yi\ 

|2 + 

1^2 - h\ 

2 ^ 1 +s 


~ h] {hihi) 


-2 


dx dy 


|mi(x) -Mi(y)p 


toi JKi, {hi |xi - 2 /i|2 + hi \X 2 - |/2p)^+" 


toxdy 


= K\ui 




\Ul 


|2 

'AHfiQ) 


|'*^l(') ^ 2 ) l_H's(0,l) 


= 


h'^j\ui{xi,X 2 ) - Ui{yi, X 2 )\‘^ dxidyi 


'0 JO 


~ h 


T WiIahukh) 


h-'-2^|Xi-|/i|l+2^ 

(*>/)> 


hi 


dxo 
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and analogously, 

(* > J)- 

Therefore, applying Theorem 15.21 fih we obtain 

||(nRTu)i||Q^^,, ~ ||(nRTu)i||g^^ < + ||divu||p^^ 

- l|willo,i^« + ^?'l“ilH=(i^«) + ^?|divu|2^^^.. (7.2) 

Similarly, applying Theorem 15.11 (i) and recalling (15.ip . we have for i > j 

||(nRTu)i||2^^,. ~ ||(nRTu)i||2^^ < 

— hihj (^ll'*^i|lo,Q + \'^^\ahi{q) I^^U^KQ)) 

— ll'^^lllgKij + Wll'AHliKij) + Wll'AH^iKij)- (7-3) 

The estimates analogous to fl7.2p and (17.3p are also valid for || (nRTu) 2 ||o,Xij with i > j. 

Combining the estimates for both components of IIrtu over all elements in T and 
then using the superposition argument as in (I5.3ip for anisotropic seminorms and the 

superadditivity property of the TT^/^-norm, we arrive at the desired result. □ 

Our construction of the operator Qh follows the technique used by Hiptmair and Schwab 
in the proof of Lemma 8.1 in |2l] but relies on stability properties of the Raviart-Thomas 
interpolation on graded meshes over individual faces of T. Given u G H^(F) fl H(div,F), 
s > 0, we consider the following mixed problem: Find (z, /) G H(div, F) x L‘l{F) such that 

(z,v) + (divv,/) = (u,v) Vv G Ho(div,F), 

(divz,^) = (divu,^) Wg e Ll{F), (7.4) 

z ■ fi = u ■ n on dF. 

Here, LiiiF) := |n G L‘^(F): (v, 1) = 0|, n is the unit outward normal vector to dF, and 
Ho(div, F) := {v G H(div, F); v fi|a^ = 0}. 

The unique solvability of (17.4p is proved by standard techniques (see [TTl Chapter II]). 
In fact, it is clear that the pair (u, 0) solves (17.4p . 

A conforming Galerkin approximation of problem (17.4p with Raviart-Thomas elements 
on the graded mesh A^(F) reads as: Find (zh, fh) G X/i(F) x Rh{F) such that 

{zh,w) + {dww,fh) = (u,v) Vv G X/j(F) n Ho(div,F), 

{dwzh,g) = {dwu,g) ^g e Rh{F), (7.5) 

z/i ■ n = Hrtu ■ h on dF. 

Here, X.h{F) denotes the restriction of X/^ onto the face F, and Rh{F) := {g G L^(F); 
g\K = const, ViC G A^(F) and {g,l) = 0}. 
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Note that the third equation in fl7.5p implies (div(u — Zh), 1) = 0. Hence, the second 
identity in fl7.5p holds for any piecewise constant function g G divX/i(F). Thus, divz/j 
is the L^(F)-projection of divu onto divX/i(F). In particular, if divu G divX/i(F) then 
div Zh = div u. 

We now prove the unique solvability of 07.51) . First, for any Qh G Rh{F) we hnd a 
function 0 G Hl{F) := {0 G (0, 1) = 0} solving the variational problem 

(V0,V0) = ((7,,0) W^eHUF). (7.6) 

Applying the standard regularity result for problem 07.61) (see, e.g., [22l p. 82]), we conclude 
that 0 G with some r G (|, 5) (here, a; < 27r denotes the maximal internal angle 

at the vertices of F), and 


||V 0 ||HqF) < \\ 4 >\\m+^iF) < Ik/ill- ( 7 - 7 ) 

Therefore, V 0 G H^(F) fl Ho(div, F), r > ^, and the interpolant HrtV 0 G X/j(F) fl 
Ho(div,F) is well dehned and stable, due to Lemma im Moreover, div (HrtV 0 ) = 

no(div V 0 ) gh- Hence, using 07 .ip and 07 . 7 p we prove the discrete inf-sup condition; 

(divv,„5f,j) (div(nRTV0),fifft) 

sup -r - r - > -f—-——^- 

v^eX;,(F)nHo(div,F) ||w||H(div,F) IIHpitV 0||H(div,F) 

> _W_ 

" c* (||V0||HqF) + ||divV0||) + ||div(nRTV0)|| 

> C\\gh\\ VghERhiF). 

This condition along with the property div(X/i(F) fl Hq (div,F)) = Rh{F) ensures ex¬ 
istence, uniqueness, and quasi-optimality of the solution {zh, fh) to 07 . 5 p (see | 11 )). In 
particular, using the quasi-optimality and recalling that z = u, / = 0, we estimate 

||u - z/j||H(div,F) < inf I|u - v,,||H(div,F) + inf \\f - gh\\ 

Vh€X/i(_F) ghGRh{F) 

{■Vh-TlKru)-h\F=0 

< ||u — nRTu||H(div,F)- ( 7 . 8 ) 

We now estimate ||u — 'Z‘h\\-R-i/ 2 (F)- One has for any e G ( 0 , |) 

|(u — Z/i, w)| 

l|n - z,,||H-i/2(n < ||u - z/,||H-i/2+e(n = sup ||—r,-. ( 7 . 9 ) 

weHl/ 2 -£(F)\{ 0 } ll’''^llHl/2-e(F) 

For a given w G we solve the following problem; Find g) G Hl{F) such that 

(V(p,V 0 ) = -(w,V 0 ) V 0 gFJ(F). ( 7 . 10 ) 
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Similarly to fl7.7l) . the regularity result for ip reads as 

ip e ^ ll/ll(_ffi/2+E(i;’))' l|w||Hi/2-e(i;’), (7-11) 

where / G is dehned by /(0) = — (w, V 0), V(/) G 

Then we set 

q ;= w + V </? G p Ho(div, F). (7.12) 

It also follows from fl7.10p that divq = divw + div V = 0. Furthermore, we have by 
(TmD-irmi) that 

||q||Hl/2-e(F) l|w||Hl/2-e(F) + \W\\h^/^-^{F) ^ 11 ^ 11 Hl/2-e (F) • (7-13) 

We now use fl7.12p and integration by parts to represent the numerator in fl7.9p as 

(u-z,„w) = (u - Z/,, q) - (u - Z/J, V <^) 

= (u - z/,, q) + (div (u - z^), ip) - ((u - z^) ■ h, 99 ) 0 ,aF- 

Hence, using (17.4p . (17.Sp and recalling that z = u, / = 0, we hnd for any qh G X/((F) (b 
Ho(div, F) and arbitrary ip^ G Rh{F) 

|(u-Z/„w)| = |(u-Z/„q-q,,) + (u-z,„q/,) 

+ (div {u-Zh),ip- iph) - ((u - Hrtu) ■ h, (^)o,aF| 

= \{u-Zh,q-qh) + (divqhjh) 

+ (div {u-Zh),ip- iph) - ((u - Hrtu) ■ h, (^)o,aF| 

< ||u-z/,|| ||q-qh|| + |(divq/„//,)| + ||div (u - z,,)|| \\ip - iph\\ 

+ ||(u — Hrtu) ■ n||F-i+£(aF) \W\\H^-^{dF)- ('^•14) 

Let denote the Raviart-Thomas interpolation operator on the ‘coarse’ quasi-uniform 
and shape-regular mesh A)('^“(F) obtained from the graded mesh A^(F) by patching to¬ 
gether long and thin elements (see Figure [3]). We also denote by Hq'^'^ the L^(F)-projector 
onto the space of piecewise constant functions on A)('^“(F). Then we set 

q, :=n^4\GX,(F)nHo(div,F) and H^J^V G R;,(F). 

By the standard properties of the Raviart-Thomas interpolation and the L^-projection 
on quasi-uniform and shape-regular meshes, we have 



div q/, 

= n;)/“divq = 

= 0, 

(7,16) 


< llq|l 

([74^ 

I|w||//i/2-6(f), 


■ q^JI 

Fl/ 2 -e(F) ^ 

(7,16) 



czm 



w^- 

^h\\ <h\\ip\\ 

773/2-6 (F) ^ 


(7,17) 
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Figure 3: An example of shape-regular quasi-uniform mesh (thicker lines) obtained by 
patching together elements of the graded mesh (thinner lines). 

To estimate || (u — IIrtu) • we recall that /g^(u — IIrtu) • h = 0 for any element 

edge Ch C dF. Therefore, we can use a standard duality argument to prove (cf. [T3l p. 259]) 

||(u-II rtu) ■ h||H-i+q9F) < ( max \eh\) ||(u - IIrtu) ■ n||R 2 (a^). 

\ehCdF / 

Then by interpolation we obtain 

II (U — IIrtu) ■ < ^ Il(u - IIrtu) ■ n|| 

H-i/2(^dF) ^ ^ ||u - IIrtuII H(div,F) 

(7.18) 

(here, we also used the continuity of the normal trace operator v i-^ v ■ hUc’ as a mapping 
H(div,F) ^ H-^/^{dF)). 

Furthermore, one has 

unj 

llv^llHi-qaF) < \\(p\\H3/2-e(^F) ~ I|w||hi/2-£(f)- (7-19) 

Now, using fl7.15p - fl7.19p in fl7.14p and recalling fl7.8p we find 

|(u - Zh,w)\ < ||u - nRTU||H(div,F) ||w||j:^l/2-e(R^). 

Using this estimate in fl7.9p we obtain 

l|u ~ ^ ||u — nRTu||H(div,F)- (7.20) 
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Now we can prove the desired result. 


Proof of Proposition 16.ll For any u G Hi(r) nH(divr, F), we define Qh^ G X/j face by 
face as Qhu\F '■= '^h for any face F C F, where z/^ is a unique (vectorial) solution to (17.51) . 
Then the commuting diagram property (16.11) follows from the second identity in (17.5p . and 
inequality (I7.2UI) yields estimate (16.21) . □ 
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